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Abstract 

This report describes the computer program GKINT, which was developed to plan, analyze 
and interpret the first High Velocity Transient Field (HVTF) ^-factor measurements on radioac- 
tive beams produced as fast fragments. The computer program and these notes were written in 
September 2004. Minor corrections and updates have been added to these notes since then. The 
experiment, NSCL experiment number 02020, Excited- state configurations in 38 S and 40 S through 
transient- field g-factor measurements on fast fragments, was performed in October 2004; results 
have been published in Physical Review Letters 96, 112503 (2006). 
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FIG. 1: Schematic view of the experimental arrangement viewed from above. Only the four SeGA 
detectors perpendicular to the magnetic field axis are shown. There are ten other SeGA detectors 
which are not shown for clarity. 

I. INTRODUCTION 

This report contains notes on the code GKINT. Much of the code was based on previous 
programs for evaluating angular correlations and kinematics in transient-field g factor mea- 
surements. These notes therefore tend to put an emphasis on the new features that have 
been introduced into GKINT. Figure ^ shows the experimental arrangement that the code 
is written to model. 

II. OUTLINE OF THE PROGRAM GKINT 

The code GKINT (G factor Kinematics at INTermediate energies) is based on previ- 
ous codes written at ANU to evaluate the quantities of interest in transient-field g factor 
measurements. The new features include (i) the use of relativistic kinematics (ii) Coulomb 
excitation calculations of cross sections and angular correlation parameters appropriate for 
intermediate energies. 
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The target can have up to 10 layers of compound composition. One layer has to be 
specified as the ferromagnetic host. In the present version of the code, the ferromagnetic 
layer is assumed to be Fe. Apart from the excitation of the 2f state of Fe, which is hard- 
wired into the code, calculations can also be performed for any other ferromagnetic host. (If 
the ferromagnet does not have Z = 26 this part of the calculation is skipped. At present the 
intermediate energy code used is not a coupled-channels code. It cannot evaluate multiple 
excitation, which might be needed for a valid calculation of the excitation of the rotational 
Gd isotopes, even in glancing collisions.) 

The main tasks performed by the code, in order, are: 

1. Read in all input data. 

2. Trace energy loss of beam through all target layers. 

3. Calculate projectile Coulex and angular correlation parameters as a function of energy- 
loss through the target at the average scattering angle. 

4. Do detailed calculations of projectile Coulex cross sections, count rates, Doppler- 
shift parameters, angular correlation parameters, transit times, ion energies/velocities 
through each layer, transient-field precessions, etc., taking into account both the finite 
solid angle of the particle counter and the energy loss of the beam in the target. 

5. Evaluate angular correlations and perturbed angular correlations, including RIV ef- 
fects, at the angles of the Segmented germanium array (SeGA). 

6. Evaluate cross sections and angular correlations for excitation of the Fe target layer. 

III. ENERGY LOSS 

A. Energy loss algorithm 

Energy loss of the beam in the target layers is evaluated using the stopping power for- 
mulation of Ziegler et al. (1985). These stopping powers appear to be quite accurate for S 
beams in Au and Fe targets. A provision is made to scale to stopping powers for each layer 
by an input parameter. 
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Because the stopping powers are a function of energy, the energy loss calculations have 
to be done by dividing the target layer into narrow strips of width Ax. To work efficiently 
with the thick foils that are used here, the energy loss algorithms were improved by using 
the Runge-Kutta method. 

To avoid confusing notation, we put dE/dx = S(E) to denote the energy dependence of 
the stopping power. If the energy of the ion as it enters the strip n is E n , then the energy 
at the next strip E n+ i is 

E n+1 =E n - AE. (1) 
We have to evaluate AE. The Runge-Kutta equations give 

AE = (AEi + 2AE 2 + 2AE 3 + AE 4 ) /6 (2) 

where 

AEi = S(E n )Ax (3) 

AE 2 = S(E n - AE 1 /2)Ax (4) 

AE 3 = S{E n - AE 2 /2)Ax (5) 

AE 4 = S(E n - AE 3 )Ax (6) 

The Runge-Kutta method gives an exact solution (but for rounding errors) when S(E) 
can be represented by a polynomial of degree not greater than four. Convergence is tested 
by repeating the calculation with double the number of steps. The energy loss converges 
within a few keV for targets several hundred mg/cm 2 thick when a few hundred steps are 
used. There can be a reduction in the accuracy if the ion almost stops, but such cases must 
be avoided in an experiment anyway. 

B. Energy loss calculations in GKINT 

The first part of the calculation is to evaluate the energy loss of the beam through the 
target layers. The beam energy at the front of each target layer is determined and the transit 
time through the layer is evaluated. The energy with which the ion leaves the target is also 
printed. 

For the subsequent integrations, each target layer is then subdivided into a number of 
strips of equal depth, and the energy of the beam at each strip is evaluated and stored. A 
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survey of the Coulomb excitation cross sections and angular correlation parameters (a 2 and 
04) as a function of depth is made, assuming that the beam scatters at the average detection 
angle, 6 = (6> max + min )/2. 

In the case of intermediate-energy fragments scattering at small angles, the difference 
between the incident and scattered beam energies is very small. However it is evaluated pre- 
cisely, and the subsequent energy loss of the scattered beam ion along its trajectory toward 
the particle counter is evaluated as a function of depth through the target, as required. 

IV. TRANSIENT-FIELD PRECESSIONS 

Transient-field precessions are evaluated for ions that traverse the ferromagnetic layer of 
the target using the routine ROTN. The primary purpose of this routine is to evaluate 



where g and r are the g factor and lifetime of the excited state. B tr is the transient-field 
strength, which varies with time as the ion slows within the ferromagnet. 

The integral over time is rewritten as an integral over energy using the non-relativistic 
relations 



where E is the energy and M is the mass of the ion. The use of the non relativistic formula 
for energy is not expected to be too unrealistic because the TF strength is small until the 
ion velocity falls to around 0.15c. 

In the routine ROTN the integral over the energy range is performed using Simpson's 
rule. Along with the precession angle, ROTN evaluates a number of other quantities that 
are often considered useful in relation to the Transient Field. 

The effective time for which the ion experienced the TF during its flight through the 
ferromagnetic foil layer must be evaluated first. From this, one may derive the average ion 
velocity and the average TF strength. 

If ti and t e are the times at which the probe ion entered into and exited from the fer- 
romagnetic foil, and r is the mean life of the probe state, the effective interaction time 
is 




(7) 




(8) 



teg = r[exp(-^/r) - exp(-t e /r)]. 



(9) 
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When t 3> U, t e , 

*eff —ti — t e — tf e , (10) 

where tf e is the time the probe ion takes to traverse the ferromagnetic layer. The average 
velocity of the ion whilst in transit through the ferromagnet is 

(v/v ) = \f\v/v Q ) exp{-t/r)dt]/t eS , (11) 

where vo = c/137 is the Bohr velocity. Again, when the mean life is much longer than the 
transit time, 

{v/vo) = L/[v t fe \, (12) 

where L is the thickness of the ferromagnetic foil. If is expressed in milliradians and the 
interaction time in ps, the average TF strength in kTesla is 

(B TF ) = O.O2O90/t cfr . (13) 

V. COULEX CALCULATIONS AND PARTICLE KINEMATICS 
A. Modifications to COULINT 

Coulex cross sections are evaluated using routines from COULINT, written by C.A. 
Bertulani, see Phys. Rev. C 68, 044609 (2003). A few modifications were needed for the 
present application. 

Firstly, the COULINT routines were modified to return the CM differential cross section 
for a specified energy and scattering angle. Both the Rutherford and Coulomb excitation 
cross sections are returned. The code was also modified to return the statistical tensor 
for the angular correlation in the projectile frame, at the specified energy and particle 
scattering angle. Since there is some difference in the terminology, and in the definitions of 
the statistical tensors, the code works with coefficients, which are the coefficients that 
occur in the familiar Legendre polynomial expansion of the angular correlation: 

W{6 y ) = 1 + a 2 P 2 (cos 7 ) + a 4 P 4 (cos 7 ) (14) 

The centre-of-mass to Lab transformations were evaluated with relativistic formulae: 

, a 1 sinflcM 

tan# Lab = — • - — 15 

72 (cos 6»cm + M P 7i/Mi7 2 ) 
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and 

— 7i72 tan fc^Lab ±1 — o tan #l ab 

Aft \ M < 7 



cos6»cm : 



(16) 



1 + 7|tan 2 6» L ab) 

where the positive sign is the only one relevant in the present work. The notation used here 
drops the prime that in many references is attached to the factors 7j. To define these gamma 
factors it is convenient to introduce p = M p /M t . As usual we have 

7 = (1 - v 2 /c 2 )- 1 ' 2 = 1 + T/m Q c 2 (17) 

where T is the beam energy in MeV per nucleon and m^c 2 = 931.494 MeV. In terms of p 
and 7 we have 

7l = 2±£ (is) 

71 (l + 2p 7 + P 2 ) 1/2 1 ' 



and 



72 = l±m (19) 

72 (l + 2p 7 + p2)i/2 ^ 



The lab cross sections were evaluated as 

(20) 



^0"Lab d& CM ^CM 



The solid angle ratio, evaluated non-relativistically in terms of the particle scattering 
angles, is 

(da/dVt) CM dQ Lah ( sin6» Lab \ 

COS(0 C M - #Lab)| (21) 



{da I 'df^Lab ^cm \sin0cM, 
The relativistic formula is 

G^Lab 72 [1 + (M p7 i/M t72 )cOs6>CM] 



(22) 



rfficM [sin 2 6 CM + 7 2 2 (cos 9 CM + M p7l /M t72 ) 2 ] 3/2 
The code uses the proper relativistic evaluation, although the relativistic formulation 
differs from the non-relativistic one by a negligible factor in many cases. 

B. GKINT description: cross sections 

Figures El and 01 show the variations in the projectile Coulomb excitation and Rutherford 
cross sections as a function of scattering angle for 38 S at 20 and 40 MeV /A, respectively, 
incident on 197 Au. The figures also show the excitation probability, Pi-*/, which is 

P^t = Jtnf )W ( 23 ) 



(da/dQ) 



Rutherford 
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FIG. 2: Cross sections as a function of scattering angle for 20 MeVA4 38 S on 197 Au. 

Note that the Rutherford cross section goes to infinity as the scattering angle approaches 
zero. The Coulex cross section remains well behaved, however, because Pi->f approaches 
zero faster than the Rutherford cross section diverges. 

GKINT performs integrals over the solid angle of the particle detector and over the depth 
of the target: 

(a-t tgt ) = 2vr / / —{9 p: E[z}) sm9 p d9 p dz, (24) 
Jo Je™ in all 

where the scattering angles and the cross section are in the lab frame, and z denotes the 
depth through the target. With the cross sections in fm 2 and target thickness in mg/cm 2 , 
the count rate R in Hz is 

R = N he ^(a-t tgt )-^ (25) 
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cross sections at 40 MeV/A 
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FIG. 3: Cross sections as a function of scattering angle for 40 MeV/^4 38 S on 197 Au. 

= (6.023 x l - 6 )iV bcam i^^l, (26) 

where iV beam is the number of beam ions/sec, Na is Avagadro's number, and A is the 
atomic mass number. The count rates for Coulomb excitation and Rutherford scattering 
as registered by the particle counter are printed for N heam = 1000. If the Rutherford cross 
section is in the region where it diverges and thus implies R > A^cam, which is clearly 
unphysical, the code puts R = Abeam and prints a warning. 
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Scattering at 3 degrees lab angle 




beam energy [MeV] 

FIG. 4: d2 angular correlation coefficient as a function of beam energy at a lab scattering angle of 
3° for 38 S on 197 Au. 

C. GKINT description: angular correlation coefficients 

The spin alignment in intermediate energy Coulomb excitation depends on both the beam 
energy and the scattering angle. As can be seen in Fig. 01 the alignment decreases as the 
beam energy decreases, eventually changing sign. Similar behavior is evident when a<i is 
plotted as a function of the scattering angle, Fig. |5J Again the alignment decreases with 
scattering angle, eventually changing sign. 

For each layer of the target, the average a k coefficients are evaluated as 

(a k ) = T^r^ I I a k (6 p ,E[z})-^-(e p ,E[z}) sm9 p d9 p dz. (27) 
\c ■ Ugt ) Jo Jo™ in ail 

The integrals are performed numerically using Simpson's rule. 

Similar integrals are performed to evaluate the average energy of excitation, the average 
transient-field precessions, etc. In each case the quantity to be averaged replaces a k in an 
expression of the same form as Eq. (j2~7l) . 
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FIG. 5: <22 angular correlation coefficient as a function of lab scattering angle for 38 S on 197 Au. 

VI. ANGULAR CORRELATION FORMALISM 

The theoretical expression for the unperturbed angular correlation after Coulomb exci- 
tation in general can be written as 

W(9 7 , 7 ) = £ J(2k + l)p kq (6 p , (f> p )F k G k Q k D%(fa, 6 7 , 0), (28) 

k,q 

where p kq {9 p ,4>p) lii the statistical tensor which defines the spin alignment of the initial 
state, and depends on the particle scattering angles (0 p ,<t> p ), F k represents the usual F- 
coefficients for the 7-ray transition, G k is the vacuum deorientation attenuation factor, 
Q k is the attenuation factor for the finite size of the 7-ray detector and -D^ o (0 7 , # 7 , 0) is 
the rotation matrix, which depends on the 7-ray detection angles (# 7 ,0 7 ). The Coulex 
calculation uses a right-handed co-ordinate frame in which the beam is along the positive 
z-axis. In the applications of interest k — 0, 2, 4. 

Since we are interested in the case where the particle detector has azimuthal symmetry, 
the statistical tensors p kq depend only on 9 P and are non-zero only when q = 0: 

PkqiOpAp) = 8q,oPko(8 P ) (29) 
The program GKINT integrates the Coulex cross sections over the opening angle of the 
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particle detector and evaluates the a k coefficients defined as 

a k = ^(2k + l)p k0 F k (30) 

In the experiments a magnetic field will be applied in the vertical direction, along the 
x-axis or = direction. We define FIELD UP as being when the external field points UP 
along the positive a>axis. The transient-field precession is given by 

A6 = g<f ) = -g^- T B tI (t)e~^ dt, (31) 

where g is the nuclear g factor, r is the meanlife of the level, and 7$ and T e are the times at 
which the ion enters into and exits from the ferromagnetic foil. (This expression corresponds 
to the precession information in the whole gamma-ray lineshape, including decays which 
occur while the nucleus is still in motion through the ferromagnet. The code also evaluates 
the precession for those nuclei that decay in vacuum after leaving the ferromagnetic layer of 
the target.) The important point to note here is that the precession angle, A9, is negative 
for a positive field (i.e. field UP) and a positive g factor. 

Now consider that the nucleus undergoes a transient-field precession of ±A9 as the ex- 
ternal field is alternatively flipped up and down. We set aside the question of whether the 
positive precession corresponds to field-up or field-down, since we don't know the sign of the 
g factor a priori. The statistical tensor is then modified by the precession about the field 
axis, according to 

p kq (±A9) = Y,Pk Q 5 Q ,oD k qQ (n/2,±A9,-n/2) (32) 
Q 

= p k0 D k q0 (n/2,±A9,0) (33) 
= p k0 D k q0 (±n/2,\A9\,0) (34) 

= pk ^wh) Y " ilAel±7T/2) (35) 

In GKINT it is convenient to work with spherical harmonics rather than the rotation 
matrices. To get the final form of the perturbed angular correlation we need to use 

0,7) = ^J^TY) Y ^^ a) (36) 

and 

Y q k *([3,a) = (-iyY k q (P,a) (37) 
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Note that, depending on the implementation in the computer routine, putting a negative (3 
value in the D matrix, and/or the spherical harmonic, may not be properly interpreted. We 
will therefore avoid doing it in the code. 

We then get the final expression for the perturbed angular correlation which GKINT 
uses: 

W(6„ 7 , ±A6) = 4vr £ j|^(-l)% fe (l A0 l ±tt/2)I* (fty, <M (38) 

The evaluation of the vacuum deorientation parameters, Gk, is discussed below. At 
present the gamma detector solid angle factors Qk are all set to unity. We will need to 
formulate a procedure to evaluate these in future (perhaps using GEANT). However the 
detectors are sufficiently far from the target that Qk = 1 is a reasonable first approximation. 



VII. LORENTZ CORRECTION FOR 7 RAYS 

Because we have azimuthal symmetry, the Lorentz correction in this case is quite straight- 
forward. The basic formulae, which have been written in many publications, are reviewed in 
the first subsection. In the second subsection the method of treating the perturbed angular 
correlation is discussed. I'm not aware of any place in the literature where this situation is 
discussed in any detail. 



A. Formulae for angular correlations 

The angle between the beam axis and the direction of emission of the photon is affected 
if the nucleus is in motion rather than at rest. The relativistic velocity addition formulae 
can be used to show that the angles of emission in the laboratory, (#i ab , 0i a b), are related to 
those in the rest frame of the nucleus, (9 nuc , (finuc), by the expressions 

cos# nuc + /3 

cos^i a b = — , 39 

1 + (3 cos 6 QUC 

and 

01ab = 0nuc, (40) 

where (3 = v/c is the velocity of the nucleus along #i ab = nuc = 0. It follows that an element 
of solid angle in the rest frame is related to an element of solid angle in the laboratory frame 
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by 

(1 + (3cos6 nuc ) 2 

To obtain the correlation in the laboratory frame, Wi a b(0iab)) from that in the rest frame of 
the nucleus, W nuc (0 n uc)5 requires both transformation of the laboratory angle to the equiv- 
alent angle in the rest frame of the nucleus using Eq. . and multiplication by the ap- 
propriate solid-angle ratio so that the 7-ray flux emitted into Atc of solid angle is conserved, 
i.e. 

/ W lah (6 lah )dQ lah = / W nnc (d nuc )dQ nnc . (42) 

This condition is satisfied if 

^lab(^lab) = ^nuc(^nuc)^^. (43) 

0.4 ilab 

B. Formulae for perturbed angular correlations 

The precession of the nucleus takes place in the frame of the nucleus (obviously!). However 
we measure the perturbed angular correlation in the Lab frame. 

The standard technique for measuring small transient-field precessions in a way that 
minimizes possible systematic errors is to form a double ratio of 7-ray intensities for each 
complementary pair of 7-ray detectors placed at angles ±# 7 with respect to the beam direc- 
tion: 



P 



N(+0 7 t) N(-fl 7 |) 
\ N(+0 7 j) N(-0 7 T) ' 1 ] 

where N(±6 I 7 f, J.) denote the numbers of counts in the photopeak of the transition of interest 

as registered in the detectors at ±0 7 for the alternative field directions (f and 1). 

The effect of the transient field is to rotate the angular correlation by a relatively small 

angle around the direction of the applied magnetic field. Thus 



dW 

W{9 1 T) = W(6^ - A0) ~ W%) - A9 



(45) 



d^ 

In the case where there are no (or negligible) Lorentz effects, the precession angle A6 is 



determined from 



A0 = ! , (46) 
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where 



and S, the logarithmic derivative of the angular correlation at the angle # 7 , is 

1 dW 



S 



W d9 



(48) 



To see how these expressions must be modified to account for Lorentz effects, it is con 
venient to note that e is formally equivalent to 



W{+64) + W{+6^) ' 
The perturbed angular correlation in the lab frame is given by 



(49) 



Hence 



W lah (6 lah T) = W nuc (0nnc T)^ = W nnc (0nnc ~ A6)^ . (50) 

W(6 nuc j) - W(6 nnc y _ 
" ab W(6 nnc [) + W(6 nuc ]) £ — 
In the case of the Lorentz boost, the expression for A6 is 

A9 = ^ . (52) 

Thus the only change in the analysis procedure is that S must be evaluated at the appropriate 
angle in the rest frame of the nucleus, which corresponds to the laboratory detection angle. 
Finally, it is worth noting that if W(8) has the form of Eq. ()14j) . then 

S(6) = -(a 2 P 21 (cos6) + a^i (cos 6))/W(6) (53) 

where P21 and P41 are associated Legendre polynomials. 



VIII. RECOIL IN VACUUM ATTENUATION FACTORS 

The Gk coefficients are evaluated as follows, based on experimental data and experience 
at ANU. 

For the experiments of interest, with ions emerging into vacuum with velocities near 
the .fT-shell electron velocity, the vacuum deorientation is dominated by the fields at the 
nucleus from the ground-state of the H-like electron configuration, i.e. Is fields. At the 

15 



lower velocities of interest 2s fields also contribute. In this description, we will begin by 
assuming that the Is fields are dominant and then introduce the 2s contribution at the end. 
The code uses both. 

Neglecting the relativistic correction (which is small for Z = 16), the field at the nucleus 
from the Is electron of an H-like ion is 

B u = 16.7Z 3 Tesla. (54) 
This free-ion field causes the nuclear spin I to precess with a Larmor frequency of 

u = g^(2I + l)Bu (55) 
The vacuum deorientation coefficients for a pure Is configuration are 

Gi * - 1 - 6 'iri£r (56) 

where r is the nuclear lifetime, and 

b k = k(k + l)/(2I + l) 2 (57) 

Thus for I = 2 we have 62 = 0.24 and 64 = 0.80. In the worst case for 1 = 2, G\ s (hardcore) = 
0.76, and although (^(hardcore) = 0.20, the k = 4 term (04) in the angular correlation is 
already small. In fact, for virtually any reasonable level lifetime (approaching a picosecond 
or more) and a g-factor value that differs even a little from zero, the hyperfine field is so big 
that G\ s « G\ s (hardcore). 

We have measured the charge-state fractions for Si and S ions with energies between 
about 90 MeV and 240 MeV emerging from Fe and Gd. (AES, Anna Wilson and Paul 
Davidson, Nuclear Instruments and Methods B 243, 265 (2006).) The results for S after Fe 
are shown in Fig. El Near v = Zvq the 15 + charge state is dominant. The other prolific 
charge states in their atomic ground states produce either smaller fields (13 + ) or no field at 
all (14 + and 16 + ). Under these conditions it is reasonable to evaluate the average vacuum 
deorientation coefficients for the experimental conditions in terms of the fraction of ions that 
carry a single electron, Qh- Thus 

G k = (1 - Q H ) + QhG\ s (58) 

After some simple algebra we get 

G k = l-Q H b k ^ T)2 (59) 
1 + (out) 2 
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FIG. 6: Charge-state fractions for S ions emerging from Fe as a function of the ion velocity. The 
highest-velocity point was obtained with beams boosted by the ANU Linac. The curves are based 
on fits to Eq. Pi . 

Thus provided that Q H < 0.5, G 2 > 0.88 and G 4 > 0.6. 

For S ions with velocities near 0.7 Zvq (100 MeV), the Li-like fields should be included: 

G k = l~ QHhj^fy, - (60) 

1 + (<jJT) z 1 + (<jJT/%y 

where Qn is the lithium-like fraction, and it is assumed that 

B ns = 16.7Z 3 /n 3 Tesla. (61) 
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TABLE I: Results of fits to Eq. Q62JI for ions with 3 or fewer electrons. 



Beam Target 


(I 


v/Zv 


range 




Parameters 




mln 


max 


(In 


t( i 




Si Fe 


14 


0.83 


1.26 


-18.6231 


26.2176 


-9.6632 




13 


0.83 


1.26 


-9.6665 


15.3430 


-6.5868 




12 


0.83 


1.26 


-3.1681 


7.0924 


-4.8990 




11 


0.83 


1.26 


-4.1570 


9.5851 


-8.0001 


S Fe: 


16 


0.67 


1.08 


-23.6376 


36.5961 


-14.9403 




15 


0.67 


1.08 


-14.9891 


26.7520 


-12.6089 




14 


0.67 


1.08 


-7.6139 


15.9538 


-9.3268 




13 


0.67 


1.08 


-3.1411 


8.4963 


-7.8757 


S Gd 


16 


0.75 


1.08 


-17.8000 


23.8431 


-8.6552 




15 


0.64 


1.08 


-12.7697 


20.9555 


-9.2347 




14 


0.64 


1.08 


-6.0770 


11.9903 


-6.6953 




13 


0.64 


1.08 


-2.2632 


5.3773 


-5.1736 


Si Gd 


14 


0.80 


1.25 


-17.9075 


23.8099 


-8.4427 




13 


0.80 


1.25 


-8.0397 


11.4628 


-4.4435 




12 


0.80 


1.25 


-5.6947 


11.2019 


-6.3201 




11 


0.80 


1.12 


0.5950 


0.1204 


-3.0000 


The experimental data for the charge 


fractions in 


Fig. IH1 were 


fitted to the form 





F(q,v/Z 1 v ) = exp(a + (Xiv/ZiVq + a 2 (v/Z 1 v ) 2 ), (62) 



where the parameters ao, ai and a 2 are determined for each charge state. The fit parameters 
for the most prolific charge states are listed in Table 

The code at present uses the fit to the charge-state results for S after Fe to interpolate 
the values of Qh and Qn and thus obtain as realistic an evaluation of Gk as possible. (It 
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will give a warning if the ion velocity is outside the region where there are experimental 
data. The extrapolations must break down because the charge-state distributions are not 
truly a gaussian function of ion velocity) 

The code uses the expression in Eq. (J6Uj) to evaluate the RIV effect. In the cases considered 
to date, lot is so large that 

(urf ^ (ujt/8) 2 _ 



1 + (cur) 2 1 + {ljt/8) 2 
and hence, to an excellent approximation, 



(63) 



G k = 1 - (Q H + Q Li )b k . (64) 

Cautionary note: We know the form of G(t) for the vacuum deorientation effect, however 
we measure 

(G k ) = (1/r) / G k (t) exp(-t/T)dt, (65) 

where t exit is the time at which the ions emerge into vacuum. The code at present assumes 
texit — in the evaluation of the (G k ), as is the case for Eq. (|56p. This gives an accurate 
result if t exit <C r. Thus the present calculations for 40 S (r ~ 20 ps) should be fine, while 
those for 38 S (r ~ 5 ps) might tend to over estimate the vacuum deorientation effect slightly. 
A more rigorous evaluation of the (Gk) value applicable in the experiments is available from 
the Monte Carlo Doppler Broadened Line Shape (DBLS) code, GKINT_MCDBLS. 



IX. SOME PEDAGOGICAL PLOTS OF THE PERTURBED CORRELATIONS 
A. TF precessions in the nuclear (projectile) frame 

While the effect of the precession is most pronounced for detectors in the horizontal plane, 
the other planes are also expected to show useful precession effects. The following figures 
[7| - IH1 illustrate how the projectile (or nuclear frame) precession varies with the phi angle of 
the 7-ray detector. 



B. RIV perturbation in the lab frame 

Fig. El shows the angular correlation in the lab frame with and without the RIV contri- 
bution. In this case the particle scattering angle is between 2.5° and 5.5°. The 38 S beam 
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FIG. 7: Perturbed angular correlations in the CM frame as observed in the horizontal, or <f> = 90°, 
plane. 
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FIG. 8: Perturbed angular correlations in the CM frame as observed in the <fi = ±48° or cj) = ±228° 
planes. 



energy is 1540 MeV. The unattenuated angular correlation coefficients are a 2 = —0.4233 and 
a 4 = -0.0990. For a H-like fraction of Q H = 0.1269 and a Li-like fraction of Q Li = 0.3108, 
G 2 = 0.8978 and G 4 = 0.6595. The Lorentz effects are evaluated using (5 = 0.089. 

Note that the RIV effect on the angular correlation is modest and that it is very difficult 
to measure the RIV attenuation because 7-ray detectors cannot be placed at either 0° or 180° 



20 



+/-100 mrad phi=37 deg 



1.5 




0.5 - 



0.0 







45 



90 

theta [degrees] 



135 



180 



FIG. 9: Perturbed angular correlations in the CM frame as observed in the (ft = ±37° or <f> = ±217° 
planes. 

to the beam in intermediate-energy Coulomb excitation measurements. The measurements 
of the charge-state fractions are therefore very important. 

X. FIGURE OF MERIT: OPTIMIZING EXPERIMENTAL CONDITIONS 

GKINT provides the quantities designated S 2 W, S 2 N and S 2 (f) 2 N evaluated at the SeGA 
detector angles as figures of merit. This section explains the rationale of these choices as 
figures of merit. 

We begin by obtaining an estimate of the experimental error on the measured precession 
for a given number of accumulated counts in a single pair of detectors. Since the precession 
angles are relatively small, we can assume that 



Thus N is the total number of counts recorded in one detector for one field direction. With 
the further approximation that 



A(+# 7 , t) ~ N{-e^ I) « A(+fl 7 , |) « JV(-0 7 , T) ~ N. 



(66) 




(67) 



it follows that 



1 
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FIG. 10: Angular correlation in the lab frame with and without the RIV contribution. 

It is usually more convenient to work with a figure of merit that scales with the number of 
counts, and hence is inversely proportional to the required beam time. Thus we can calculate 
quantities proportional to S 2 N as the figure of merit and note that we must maximize it for 
optimum sensitivity. 

GKINT calculates both S 2 W and S 2 Wa, denoting the latter as S 2 N. Generally S 2 N = 
S 2 Wa will be the quantity of interest. (S 2 W is retained in the output to give an indication 
of the angular correlation contribution, separate from the cross section.) These figures of 
merit are independent of the magnitude of A9. Clearly, in a real experiment the goal is to 
minimize a^ e / 'AO. As another figure of merit GKINT also evaluates 

S 2 <p 2 Wae^ T ^ /T \ (68) 

which is proportional to (a&o/ AO) 2 for the fraction of excited nuclei that experience the full 
transient-field precession and decay in vacuum beyond the target. Here T cxit is the average 
time at which beam ions exit from the target. 

The calculations are given with and without the effect of the RIV attenuation. Since the 
RIV effect can not be avoided, the unattenuated angular correlations are given mainly as a 
point of reference. 
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TABLE II: Additional details of the angular correlation in Fig ^1] for the case without vacuum 
deorientation. 



posn 


e 


6 

T 


Chi finnr / Cfci ilnh 




S(0) 


S 2 W 


1 


24.0 


180.0 


1.1763 


0.8013 


0.0000 


0.000 


2 


40.0 


131.0 


1.1433 


1.0459 


0.5937 


0.369 


3 


90.0 


112.0 


0.9920 


1.1630 


-0.0381 


0.002 


4 


147.0 


143.0 


0.8584 


0.6329 


-0.6467 


0.265 


5 


139.0 


228.0 


0.8705 


0.7315 


0.6896 


0.348 


6 


60.0 


241.0 


1.0870 


1.2077 


-0.2589 


0.081 


7 


151.0 


90.0 


0.8533 


0.5861 


-1.1148 


0.728 


8 


151.0 


270.0 


0.8533 


0.5861 


1.1148 


0.728 


9 


29.0 


91.0 


1.1675 


0.8816 


1.0543 


0.980 


14 


29.0 


271.0 


1.1675 


0.8816 


-1.0543 


0.980 


17 


24.0 


1.3 


1.1763 


0.8013 


0.0254 


0.001 


18 


60.0 


61.0 


1.0870 


1.2077 


0.2589 


0.081 


19 


139.0 


46.0 


0.8705 


0.7315 


-0.6675 


0.326 


20 


147.0 


323.0 


0.8584 


0.6329 


0.6467 


0.265 


21 


90.0 


292.0 


0.9920 


1.1630 


0.0381 


0.002 


22 


40.0 


310.0 


1.1433 


1.0459 


-0.6027 


0.380 



XI. EFFECT OF A MASK ON THE PARTICLE DETECTOR 

Table IIVI summarizes some calculations made to assess the use of different masks on the 
particle detector. The main benefit of the mask is to reduce the Rutherford rate. The small 
reduction in Coulex count rate is compensated for by an increase in the anisotropy of the 
angular correlation. 



23 



TABLE III: Additional details of the RIV attenuated angular correlation in Fig llUl 



posn 9 <j> dn nuc /dQ lSbh W(0) uh S(8) S 2 W 



1 


24.0 


180.0 


1.1763 


0.8450 


0.0000 


0.000 


2 


40.0 


131.0 


1.1433 


1.0454 


0.5165 


0.279 


3 


90.0 


112.0 


0.9920 


1.1536 


-0.0468 


0.003 


4 


147.0 


143.0 


0.8584 


0.6562 


-0.5152 


0.174 


5 


139.0 


228.0 


0.8705 


0.7404 


0.5773 


0.247 


6 


60.0 


241.0 


1.0870 


1.1911 


-0.2669 


0.085 


7 


151.0 


90.0 


0.8533 


0.6169 


-0.8662 


0.463 


8 


151.0 


270.0 


0.8533 


0.6169 


0.8662 


0.463 


9 


29.0 


91.0 


1.1675 


0.9096 


0.8472 


0.653 


14 


29.0 


271.0 


1.1675 


0.9096 


-0.8472 


0.653 


17 


24.0 


1.3 


1.1763 


0.8450 


0.0196 


0.000 


18 


60.0 


61.0 


1.0870 


1.1911 


0.2669 


0.085 


19 


139.0 


46.0 


0.8705 


0.7404 


-0.5588 


0.231 


20 


147.0 


323.0 


0.8584 


0.6562 


0.5152 


0.174 


21 


90.0 


292.0 


0.9920 


1.1536 


0.0468 


0.003 


22 


40.0 


310.0 


1.1433 


1.0454 


-0.5242 


0.287 
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TABLE IV: Effect of a mask on the particle detector. The calculations are for a standard S 
configuration with .Ebeam = 1540 MeV. The target is 355 mg/cm 2 Au followed by 100 mg/cm 2 . 
The stopping powers are scaled by a factor of 1.02. 

6f n Count rate/1000 beam ions (Hz) S 2 Wa 

Coulex Rutherford 



#£ iax = 5° 

0.0596 1000 762 

1 0.0598 290 778 

2 0.0523 63 856 
2.5 0.0466 36 869 

3 0.0395 21 836 

= g go 

0.0720 1000 1119 

1 0.0713 292 1138 

2 0.0647 66 1236 
2.5 0.0589 38 1260 

3 0.0518 24 1237 
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